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Abstract 

Near threshold neutral pion electroproduction on the deuteron is studied in the framework of baryon chiral 
perturbation theory at next-to-leading order in the chiral expansion. We develop the multipole decom- 
position for pion production off spin-1 particles appropriate for the threshold region. The existing data 
at photon virtuality k 2 — —0.1 GeV 2 can be described satisfactorily. Furthermore, the prediction for the 
S-wave multipole Ed at the photon point is in good agreement with the data. 
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1 Introduction 



Pion photo- and electroproduction off single nuclcons in the threshold region can be considered one of the best 
testing grounds for our understanding of the chiral pion-nucleon dynamics resulting from the symmetry structure 
of QCD (for a recent status report, see e.g. Q). In the absence of neutron targets, it is mandatory to consider 
pion production off light nuclei which also leads to the consideration of interesting aspects related to few-nucleon 
dynamics. In this paper, we consider pion electroproduction on the deuteron above threshold extending our 
previous work ||. This is mandated by the following developments: First, the threshold results obtained in 
Q can not be directly compared to the data. Furthermore, the important single scattering contribution was 
not calculated to fourth order (which is mandatory to describe the elementary process with sufficient accuracy) 
but simply shifted to its value at the photon point. Such a procedure is only well controlled for the transverse 
multipole. Second, coherent neutral pion production off deuterium at photon virtuality k 2 — — O.lGeV 2 has 
been measured and analyzed at MAMI ||], and these data show a significant discrepancy in the dominant 
longitudinal cross section (S-wave multipole) from the prediction of 0. In addition, new measurements of 
pion electroproduction off the proton at MAMI Q at the lower photon virtuality of k 2 = —0.05 GeV 2 have 
led to intriguing results that can neither by explained in chiral perturbation theory nor with any sophisticated 
model (note the very unusual values for certain P-waves given in that paper). Here, we want to improve the 
calculation for coherent pion production off the deuteron in two ways. First, in the single scattering contribution 
we include the full fourth order result for the transverse and longitudinal S-waves, with its parameters fixed 
from recent data on neutral pion photoproduction and the older NIKHEF || and MAMI measurements 
for electroproduction off the proton at k 2 = —0.1 GeV 2 . Note that it was recently shown that in the case 
of neutral pion photoproduction the fourth order corrections to the P-wave multipoles are fairly small jjj. 
A similar analysis for electroproduction is not yet available, it is, however, conceivable that similar trends 
will persist in that case. Second, we calculate above threshold, which leads to a considerable complication 
in terms of the multipole expansion. While this formalism has already been developed in js| [9J, we present 
here a new form particularly suited for the threshold region and that most closely resembles the single nucleon 
multipole expansion. We restrict ourselves to S- and P-waves and evaluate the three-body corrections (or 
meson-exchange currents) to third order in the chiral expansion. We include, however, the pion mass difference 
which is formally of higher order but constitutes the dominant isospin breaking effect. The single scattering part 
for the proton has been fixed before, and can be reliably estimated for the neutron using resonance saturation 
at the photon point. However, at finite photon virtuality the situation is less clear and we do not want to rely 
on the resonance saturation hypothesis. We therefore perform two types of fits. In the minimal fit we employ 
resonance saturation for a dimension four LEC (there are in principle two LECs but their sum is constrained 
by a low-energy theorem |l(|) and use the longitudinal deuteron multipole Ld as extracted from the MAMI 
data to pin down one parameter related to a particular dimension five operator p0| . In a second scenario, 
we do not use resonance saturation and thus have two free parameters related to the polynomial part of Lq, 
which we determine from a best fit to the measured total cross sections at low pion excess energies. This still 
leaves sufficient predictive power since we can compare directly with the measured differential cross sections or 
the extracted S-wave cross section a^d- For doing that, one has to select a deuteron wave function. In we 
employed the hybrid approach using various high precision wave functions together with the chirally expanded 
interaction kernel. Here, we also improve on that aspect using recently obtained precise effective field theory 



wave functions that are consistent with the power counting of the kernel JO], [lj, [T3|,^ We will demonstrate 
that this improved calculation is in fair agreement with the MAMI deuteron data at k 2 = —0.1 GeV 2 ||, thus 
solving one apparent discrepancy and deepening the mystery surrounding the data of Ref . [|| . Needless to say 
that a separate investigation of this second puzzle is urgently called for but should not be a topic of the present 
paper. 

This paper is organized as follows. In Section @ we discuss the multipole decomposition for neutral pion 
electroproduction off a spin-1 target. In Sectionl3 we briefly review the effective Lagrangian underlying the 
calculation and the standard power counting formulas. In Section ^ the calculation of the various contributions 
to the transition current (single scattering and three-body terms) is outlined. Section || contains the results and 
discussions thereof. A brief summary and outlook is given in Section [| The appendices include our conventions 
and give many more details on the calculations. 



# 2 We have also performed calculations using precise phenomenological wave functions as a check. None of the results shown later 
depend on the choice of wave function. 
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2 Multipole decomposition 



The main part of this section is concerned with the multipole decomposition for the process j*d — + Tr°d. 
To develop this, we heavily rely on the work of Arenhovel || ^ for the classification of the operator basis, 
construction of invariant amplitudes and the calculation of observables. However, in his work the main emphasis 
was put on the helicity basis. A formal proof of the equivalence between the multipole expansion used here and 
the one of Arenhovel is given in appendix [A|. We also summarize some basic formulae to calculate observables 
from the multipoles. 

The invariant matrix element for the process j*(k) + d(pd) — > n°(q) + d(p' d ), where 7* denotes the virtual 
photon with virtuality k 2 < 0, d the deuteron and 7r° the neutral pion with four-momentum = (u>,q), can 
be expressed in terms of 13 invariant functions, 
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(2.1) 



where the 13 operators 0\ are expressed in terms of combinations of the direction of the photon three- 
momentum k, the photon polarization vector e x , the direction of the pion three-momentum q and the deuteron 
spin vector S. Here, A denotes the helicity of the in-coming photon, with A = 0,±1. This non-relativistic form 
is most appropriate for near threshold production. The explicit form of the O x , first written down in is 
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the symmetric traceless tensor of second rank in standard notation. The first nine of these operators are 
transverse, whereas the other four are longitudinal. Note also that we work in the Coulomb gauge Sq = 0. 
Furthermore, the Fi are functions of three kinematical variables, we chose here to work with the pion energy u), 
the photon virtuality and the scattering angle, Fi = Fi(u>, fc 2 , z). More precisely, 9 is the scattering angle in the 
ir°d center-of-mass system with z — cos 9. In what follows, we will however not display these arguments. For 
later use, we define the vector F via 

F = (F 1 ,F 2 ,F 3 ,...,F 13 ) . (2.5) 

Any given tree or loop graph can now be expanded in this basis, and all observables can be expressed as functions 
of the Fi . Explicit expressions can be found in || . 

However, for the analysis of the data and the direct comparison with theoretical predictions, it is advantageous to 
use a multipole decompositon similar to the standard case of pion production off a single nucleon. The pertinent 
method to do that has been outlined a long time ago in juj, and we use that formalism here to develop the 
multipole decomposition for our case. As shown in Fig. [l], a photon with helicity A and multipolarity L produces 
the neutral pion. In the final-state ir°d system, the pion has relative orbital angular momentum L v , which 
couples with the deuteron spin to the total angular momentum J. In the ir°d frame, the invariant matrix 
element for the unpolarized case considered here can be written as 



M X = \m')(qm'\a x J^(uj,k 2 )\k\m)(Xm 



3 



L 



71 



L 




Figure 1: Graphical representation of the angular mo- 
menta involved in pion electroproduction. The in-coming 
photon (wiggly line) has helicity A and multipolarity L. The 
final ir°d state is characterized by the total angular momen- 
tum J and the pion (dashed line) angular momentum L^. 
The double line denotes the deuteron. 
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where = — (eg /fen) is the transverse polarisation vector and is the photon vector current impinging 
on the deuteron and we have summed over the initial and final- state deuteron magnetic quantum numbers. 
This matrix element can be expressed in terms of electric E, magnetic M and longitudinal L multipoles ag^ 
(note that we do not give the explicit dependence on the azimuthal angle ip here since we will only consider the 
transverse and the longitudinal cross section in what follows, see also the discussion below) 
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The first factor is proportional to the angular momentum eigenfunction in the final state (i.e. the appropriate 
spherical harmonics), the second term collects the pertinent Clebsch-Gordan coupling coefficients and the third 
term contains the dynamical information in terms of the multipoles. These multipoles depend on the pion energy 
and the photon virtuality. They are characterized by three labels. The superscript L refers to the multipolarity 
while the lower indices L„ and J denote the orbital angular momentum and the total angular momentum of 
the final pion-deuteron system. We furthermore use A = \J2A + 1 for angular momentum eigenvalues. Note 
also that for real photons with k 2 = (photoproduction), the photon has no l ongitudinal components and thus 
there is no coupling to the longitudinal multipoles. One can invert Eq. (2/7) and project out the multipoles, 
this gives 
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where due to parity, the sum L+L w has to be odd for the electric and the longitudinal multipoles and even for the 
magnetic ones. Since now for a given orbital angular momentum we have the conditions |L W — 1\ < J < + 1 

* 3 Note the symbol L is used for the multipolarity and the longitudinal multipoles. However, no confusion can arise since it is 
always obvious from the context what is meant. 
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and | J — 1| < L < J +1, parity allows for four different electric, four longitudinal and five magnetic multipoles, 
which we collect in the nine-component transverse vector Tl^, 
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and the four-component longitudinal vector Ll , 



These together define the multipole vector Ml„i 
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which has 13 components. It is straightforward albeit somewhat tedious to work out the transformation matrices 
between the multipole basis and the one spanned by the invariant functions F^ . The multipoles can be obtained 
from the Fi by a 13x13 block-diagonal matrix, such that 
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where Dl„ is a 9 x 9 and El„ a 4 x 4 matrix. The explicit representation of these matrices in terms of Legendre 
polynomials is given in appendix [b|. Similarly, the inverse transformation is given in terms of a 9 x 9 matrix, 
called Gl„ and a 4 x 4 matrix, denoted Hl , as 
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The explicit form of these matrices is also given in appendix M As a non-trivial check we have shown that the 



product of the two 13 x 13 matrices in Eqs. (2.13,2.14) is indeed the unit matrix. Note that for the electric and 
the magnetic multipoles L has to be larger or equal to one (since |A| = 1) and that — or L — implies 
J = 1. This reduces the number of allowed multipoles for L n < 2 and leads to the lowest permissible value of 
L„ for the various multipoles given in Table 0. 
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Table 1: Lowest permissible value of L x for the various multipoles. 



The unpolarized differential cross section for neutral pion electroproduction from a spin-1 target decomposes 
into four structure functions. However, so far no data are available for the small transverse-longitudinal and 
transverse-transverse interference structure functions and we thus will also not consider any azimuthal depen- 
dence here (i.e. angular dependence between the scattering and the production plane). Then, the differential 
cross section decomposes into a transverse and a longitudinal part, the latter being multiplied by el, with el 
the longitudinal degree of virtual photon polarization which is related to the transverse one by 

k 2 

£ l = -T2 £ > ( 2 - 15 ) 

Kr, 
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with the photon energy and momentum taken in the photon-deuteron center-of-mass system. The explicit form 
of the transverse cross section reads 
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Note that there are no electric times magnetic multipole interference terms due to the selection rules given 
above. Similarly, the longitudinal cross section is given entirely in terms of the longitudinal multipoles 



^ = jfj \ E / dQ K'fi, m (9f 
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At threshold, only the three multipoles -Bqij M?i an d An contribute and these define the transverse and 
longitudinal S-wave cross section aod as used in M , 

a d = \E d \ 2 +e L \L d \ 2 , (2.18) 

with | £7 d | 2 = |-Eqi| 2 + lM?il 2 an d \Ld\ 2 = l^oil 2 - This concludes the formalism needed in this paper. 



3 Effective field theory 

In this section, we briefly discuss the effective chiral Lagrangian underlying our calculations and the correspond- 
ing power counting. For previous related work on pion photoproduction off nuclei see Jl5| , |l6|] . 

At low energies, the relevant degrees of freedom are hadrons, in particular the Goldstone bosons linked to the 
spontaneous symmetry violation. We consider here the two flavor case and thus deal with the triplet of pions, 
collected in the matrix U (x). It is straightforward to build an effective Lagrangian to describe their interactions, 
called /Ctttt- This Lagrangian admits a dual expansion in small (external) momenta and quark (meson) masses 
as detailed below. Matter fields such as nucleons can also be included in the effective field theory based on 
the familiar notions of non-linearly realized chiral symmetry. These pertinent effective Lagrangian splits into 
two parts, CnN and Cnn, with the first (second) one consisting of terms with exactly one (two) nucleon(s) 
in the initial and the final state. Terms with more nucleon fields are of no relevance to our calculation. The 
pertinent contributions to neutral pion photoproduction at threshold are organized according to the standard 
power counting rules, which for a generic matrix clement involving the interaction of any number of pions and 
nucleons can then be written in the form 

M = q v T(q/y), (3.1) 

where /i is a renormalization scale, and v is a counting index, i.e. the chiral dimension of any Feynman graph. 
v is, of course, intimately connected to the chiral dimension di which orders the various terms in the underlying 
effective Lagrangian (for details, see [f?|). For processes with the same number of nucleon lines in the initial 
and final state (A), one finds 

v = 4 - A-2C + 2L + E^ A i 

i 

A t ee d t + n,/2 - 2. (3.2) 

where L is the number of loops, V% is the number of vertices of type i, is the number of derivatives or powers 
of M n which contribute to an interaction of type i with rii nucleon fields, and C is the number of separately 
connected pieces. This formula is important because chiral symmetry places a lower bound: Aj > 0. Hence 
the leading irreducible graphs are tree graphs (L = 0) with the maximum number C of separately connected 
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G 



pieces, constructed from vertices with Aj = 0. In the presence of an electromagnetic field, this formula is 
slightly modified. Photons couple via the electromagnetic field strength tensor and by minimal substitution. 
This has the simple effect of modifying the lower bound on A; to A; > -1, and of introducing an expansion in 
the electromagnetic coupling, e. Throughout, we work to first order in e, with one exception to be discussed 
below. For more details on the counting, we refer to Eg ], In what follows, we will work within the one- 
loop approximation to order q 3 (notice that we refer here to the chiral dimension used to organize the various 
terms in the calculation of the single-nucleon photoproduction amplitudes), with the exception of the S-wave 
contribution to the elementary process j*N — ► 7r°iV (as discussed in the introduction). In terms of the counting 
index v, we include all terms with v = 4 — 3A = —2 and v = 5 — 3A = — 1. Consequently, the effective 
Lagrangian consists of the following pieces: 

r _ r{2) , r (l) , ~(3) r , ~(4) i , AO) , r (2) /„ o\ 

•'-eft — 7T7T ' '"ttN ' '"itN ' '"itNI'^kNI ' *"NN ' '"NN ' \°-°) 

where the index (i) gives the chiral dimension di (number of derivative and/or meson mass insertions). The 
form of C^tt +£^. 1 jy is standard. The terms from C^+C^Jj contributing to the single-nucleon electroproduction 



amplitudes are given in Ref.EjJ. Note that the square brackets in Eq. ( |3.3[ ) indicate that such fourth order 



terms are only taken for the S-wave single nucleon production amplitudes^ 4 . The effective Lagrangian can also 
be used to generate deuteron wave functions of sufficient precision, as done in Refs. Jl2] , [i"3| based on a modified 
Weinberg power counting. We use the NLO and the node-less NNLO* wave functions from 13 for the allowed 
cut-off range A = 500 . . . 600 MeV, where the cut-off stems from the regulator function in the Lippmann- 
Schwinger equation used to generate the bound and the scattering states (for a more detailed discussion giving 
also the explicit form of C^ff , see e.g. Ref. |ll]]). As a check, we have also made use of the hybrid approach 
of |l!| , sewing precise phenomenological wave functions to the chirally expanded kernel. None of the results 
discussed later depend on the choice of wave function and we therefore will only present numbers for the chiral 
EFT wave functions. After these general remarks, let us now turn to the actual calculations. 



4 Anatomy of the calculation 

In this section, we outline how the various contributions to the multipoles and the observables are calculated. 
First, we briefly discuss the separation of the transition matrix elements into two- and three-body terms (or, 
in nuclear physics language, impulse and meson-exchange terms). Then, these two types of contributions are 
discussed separately, in particular we stress the differences to the threshold calculation of [0 . Many details are 
relegated to the appendices. 

4.1 General remarks 

Consider first a generic diagram for neutral pion electroproduction off the deuteron as shown in Fig. |^. The 
interaction kernel decomposes into two distinct and different pieces. First, the virtual photon can produce the 
pion on either the proton of the neutron, with the other nucleon acting as a mere spectator. This is called the 
single scattering contribution (ss), compare Fig. |^. It is important to note that to properly account for this one 
has to transform from the photon-nucleon to the photon-deuteron center-of-mass system as discussed below. 
Second, all other terms in the interaction kernel involve both nucleons, comprising the so-called three-body 
(tb) interactions (see again Fig. |^). To be specific, the transition matrix element for pion production of the 
deuteron has the form 

M x = Y, \0(p d 'qK^\p d k)(m s \ (4.1) 

m s ,m' s 

in terms of the three-momenta of the in-coming (out-going) deuteron, pd and p^' , respectively. Here, J% 
denotes the vector current impinging on the deuteron. To unravel the underlying structure of the deuteron, one 
expresses this matrix element in terms of the two-nucleon current, Jnn- This current then separates into the 
two terms just discussed, 

{p d 'q\a^\p d k) = ^^ 4E ^ E > 2 {Pd ' ^ J " NN Wd % ] 

* 4 In fact, one also has to account for one particular dimension five operator as explained in |10[ , see also Sect. F'd 
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tb 



Figure 2: Decomposition of the full interaction kernel (as shown in 
leftmost diagram by the shaded circle) into the single scattering (ss) and 
the three-body (tb) contribution. The triangle symbolizes the deuteron 
wave function. 



where Ei (E^) denotes the energy of nucleon i (i = 1, 2) in the initial (final) state. In what follows, we will work in 
the approximation that pion is produced in S- or P-waves only, that is we allow for — 0, 1. This means that we 
have to consider three S-wave (E^ v Lqd Mn) an d eight P-wave multipoles {E\ l ,Ef 2 , M\ , M-jj, M^ 2 , Mf 2 , L° v 
Lf 1; L\ 2 ), compare Table [jj Only if one allows for higher partial waves, all thirteen different structures discussed 
in Section ^ will contribute. We will now discuss the single scattering and the three-body contributions in more 
detail. 



4.2 Single scattering contribution 

In the previous paragraph, we have introduced the single scattering contribution to the two-nucleon current. 
It has the following generic form 



i(2m(27r) 3 5(pi'-pi) 



A jtt PlM + a X r n„ 

(J, fJi 



(1 



(4.3) 



with pi — pi ' = p — p' — k/2 + q/2 in terms of the in-coming (out-going) relative nucleon cms momenta p(p'), 
the pion and the photon momenta and similarly for pi — pi' ■ The isoscalar current is then expressed in terms 
of the conventional CGLN |20t] amplitudes, 
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(4.4) 



where the stared quantities refer to the center-of-mass system of nucleon i, the M£ s denote the single nucleon 
multipoles, the explicit form of the operators Of s is given in appendix |^ and the corresponding transformation 
matrix for the transverse and longitudinal multipoles is standard pl| . Of course, one has to transform these 
expressions from the pion-nucleon to the pion-deuteron center-of-mass system. This is described in some detail 
in appendix When sandwiched between the deuteron wave functions, one has to deal with integrals of the 
type 



d 3 p ( 



p -2 + 2 



O x + 2 -p) 4>{p) 



(4.5) 



where 4>(p) denotes the momentum space deuteron wave function, and 0%, O2 are arbitrary spin structures. To 
work in coordinate space, one has to Fourier transform these expressions, using the basic integrals 



d 3 r</» t (f)O 1 0(f)e- l (' ; ~ 9 "')- r "' /2 , 
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dV* Ip - - + f j 2 , iPi d> [p) = J d 3 r 0t(f) G 2il i 9 ri 0(f) e-<^«)' ? / 2 , (4.6) 

and further decomposing the deuteron co-ordinate space wave function 4>{f) into its radial S- and D-wave 
components, 

ir \ 1 I u ( r ) flwfr) „ , A .\ , . 



4tt \ r V 8 r 

Here, $12 (f) is the usual second order tensor operator. We also remark that in [B we had expressed the overlap 
integrals in a factorized form, more specifically, the transverse and longitudinal deuteron multipoles could be 
written as products of the single nucleon S-wave multipoles and a set of deuteron form factors. Such a procedure 
becomes very complicated above threshold and is not transparent, therefore we do not follow such a path here 
in detail (although we have performed some calculations in that framework to have a further check on the 
numerics). Still, it is important to stress that the single scattering contribution is strongly suppressed with 



increasing photon virtuality because of the decreasing overlap integrals in Eq. (4.5) 



Again, we work in the S- and P-wave approximation. The corresponding single nucleon multipoles are subject 
to the chiral expansion. We work to first non-trivial loop order, i.e. to third order, with the exception of 
the proton and neutron S-wave multipoles -Bg-f an< ^ • We take the form given in |Tcj which includes all 
fourth order terms and one particular fifth order term necessary to separate cleanly the longitudinal from the 
transverse piece. We refer to that paper for the explicit expressions of the single nucleon transition current. 
Here, we only spell out the generic form for the S-waves, 

g _ ^Born + £q 3 -loop + 5 q 4 -loop + get ^ ^ 

with S = Eq^ 1 or £q+. At fourth order one has two local operators ~ k 2 with the LECs a|, a\ (I — p,n). 
However, a particular low-energy theorem (LET) strongly correlates these two LECs, in the soft-pion limit one 
has a^ + a^ — 0. To break this correlation that is not observed in the proton data, one has to include a correction 
to the LET of the form Lq^I 1 = —eM 2 k 2 a\ which is formally of fifth chiral order (in the effective Lagrangian). 
As noted before, with that we are able to describe the proton data of ||, |) for j*p — > ir°p at k 2 — —0.1 GeV 2 
but not the more recent data of at half the photon virtuality (as discussed in the introduction). For the LECs 
related to the neutron amplitude, we follow two strategies. First, we fix — a J from resonance saturation as 
detailed in Ref. (lO) and determine a% from a fit to the empirical threshold amplitude Ld of Ref . || . Second, we 
leave both and as free parameters thus relaxing the constraint due to the LET and fit to the threshold 
total cross sections of Ref. ||. In what follows, we will call these two procedures fit 1 and fit 2, respectively. The 
scaling procedure performed in Q was done too simplistically for the longitudinal S-wave multipole leading to 
the too large S-wave cross section. 

There is one additional point that deserves particular discussion. With increasing photon excess energy (that is 
the energy normalized to the threshold energy), one should observe two unitary cusps due to the opening of the 
n + nn and the ir~pp channel for neutral pion production off the proton and the neutron, respectively. Indeed, if 
one does not boost the energy dependent multipoles as suggested in p2|] , two cusps are visible as shown in Fig. || 
for the real part of E\ together with the corresponding growth of the imaginary part^. The effect of applying 



the full boost correction, i.e. also changing the arguments of the multipole amplitudes, turns out to be small, 
as also shown in Fig. ^. The real part is shifted by a few percent and the imaginary part is almost unaffected 
for the pion energies considered here. Note that in all non-analytic terms, like e.g. the energy dependence of 
the S-waves ~ y/l — uj 2 /u> 2 (and similar terms in the P-waves), we always use the physical values for the nnir + 
and ppn~ thresholds given by uj c . This is consistent with the chiral expansion and the analytic structure of the 
amplitudes. 

4.3 Three-body contribution 

We now turn to the three-body contribution. Above threshold, we have 8 diagrams contributing at third order, 
see Fig. ^. Note that in contrast to earlier works, we differentiate between the charged and neutral pion masses 

# 5 Note that we work with AI^o = 134.97 MeV and M n+ = 140.13 MeV or = 142.53 MeV, to account for the neutron-proton 

mass difference in the rescattering diagrams. A detailed discussion of this point is given in Ref. |23| . 
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Figure 3: Effect of the Lorentz boost on the arguments of the multipolcs. 
Shown are the real (upper two lines) and the imaginary part (lower 
two lines) of the single scattering contribution to the electric multipolc 
for varying pion excess energy AW at fixed photon virtuality k 2 = 
—0.1 GeV 2 . Solid (dot-dashed) lines: with (without) boost. 



also in these graphs, although this is formally an effect of higher order. We remark that one can combine these 
various contributions into two distinct classes with one and two pion propagators, respectively, the so-called 
rescattering and pion-in-flight diagrams. 



a) 



b) 



c) 



d) 



c) 



f) 



g) 



h) 



Figure 4: Three-body interactions which contribute to neutral pion 
electroproduction at threshold to order q 3 (in the Coulomb gauge). Here, 
solid, dashed and wiggly lines denote nucleons, pions and photons, in 
order. 

The explicit analytical expressions for the corresponding matrix elements in momentum space are: 
Diagrams a) + b) + c) + d) + c) (rescattering type) 



■Ma+b+c+d+e 



F 3 



£■ S 



2 1 - ^ 

9 a — 1-1 -Qo 
9o 



- g A — S ■ q e ■ q ) 

qo J q 2 + 5 2 



[fl • T 2 - T 3 T 3 } 



(4.9) 



Here, q' = {q' 0l q ) with q' = q + 0(l/m) is the four-momentum of the exchanged pion and 



5*=Ml + -ql-ie. 



(4.10) 
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Diagrams f) + g) + h) (pion- in- flight type) 



M f+g+h = -2- Fr e.( q +q)S-q r2 + M 2 +)((f2 + a 2) [n-rz-T^] . (4.11) 

Here, we use the following convention. The intermediate pion has momentum q after emission from the left 
nucleon line and before the interaction with the photon. After that, the momentum is q and the pion is 
absorbed on the right nucleon. Note that for both classes of diagrams, the factor 2 in front takes care of the 
interchange 1^2. We remark that in contrast to previous work ^ [T(| we differentiate between the charged 
and neutral pion masses for the exchanged meson. While that is formally an effect of higher order, we still 
consider it here because it was already shown to be the dominant isospin breaking effect in the investigation 
of pion photoproduction off nucleons, first discussed in the context of chiral perturbation theory in [j24| . In 
essence, we have calculated all three-body graphs for the two different values of to c corresponding to the opening 
of the ppn~ and the nmr + channels and performed the necessary average. Since these operators are sandwiched 
between the deuteron wave functions and one has to integrate over all momenta, one picks up an imaginary 
part from the intermediate NNtt state, i.e. the corresponding propagators have to be split into a real (principal 
value) and an imaginary part. It appears when the momentum of the exchanged pion is equal to the charged 
pion mass, that is at an excess energy of 



AW C =lu c + ^m 2 d + ujI - M 2 - Wo = 5.3 (7.9) MeV , (4.12) 

for the ppir~ (nmr + ) intermediate state. Here, Wo ~ rrid + M T o is the threshold energy. These poles will reveal 
itself as a cusp-like structure in the corresponding multipoles (if one considers the three-body terms separately, 
see Section ||). As before, we have restricted the calculation to relative S- and P-waves. We have also performed 
calculations without this truncation, which gives a relative measure of the contribution from D- (and higher 
partial) waves. More specifically, consider a typical coordinate space integral, 

oo 

dn r e-<^/ 2 y^47rj (ar) = 4tt ^ (-sgn( 2 ;)) i (2L + 1) ^(^(c) F L (g ■ fc) , (4.13) 

with a = \kx + <f/2| r, b = kr\x\ and c = qr/2. Note that the coefficient a in the Bessel function jo(o-r) contains 
the explicit angular dependence. In case of the S- and P-wave approximation, one operates with the projection 
matrix on the full sum and the series in Eq. ( [4.13 ) truncates after the first few terms. We have found that these 



differences are very small, justifying a posteriori the assumption of only retaining the S- and P-waves. 



5 Results and discussion 

In this section, we display the results for the multipoles, differential and total cross sections and the S-wave 
cross section aod for the two fit strategies. We have performed calculations with the chiral EFT wave functions 
at NLO and NNLO* Jlf , [[J) for cut-offs in the range of 500 to 600 MeV. Since the results for the observables 



are very similar for all these various wave functions, we only show these for the NNLO* wave function with 
A = 600 MeV. The fitted LECs vary mildly for the various wave functions as shown in Table ||. We note that 
all these LECs are of natural size (note that the value for a§ appears unnaturally large due to the particular 
definition of this LEC as used in Ref . JlOtl , see also the discussion in that paper). We remark that the results 
using the Bonn wave function as employed in Ref. Q are fully consistent with the ones based on the chiral EFT 
wave functions and we thus refrain from displaying these numbers here. 



The real parts of the various S- and P-wave multipoles are displayed in Figures |5j-|§| by the solid lines for fit 2 
using the NNLO* wave functions with A = 600 MeV. The corresponding single scattering contribution is also 
shown (dashed lines) ^\ As found in previous calculations, the three-body effects are sizeable, especially in 



the S-waves. In contrast to previous attempts using meson-exchange models this does not pose a problem for 
extracting the single scattering contribution because one can systematically calculate the higher order corrections 

# 6 We refrain from showing these multipoles using the other wave functions for fit 2 because they come out very similar. There 
are some differences in the multipoles for fit 1, as will be discussed for the S-waves later on. 
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Figure 5: S-wave multipoles as a function of the pion excess energy AW 
for photon virtuality k 2 = —0.1 GcV 2 . The solid (dashed) line shows the 
total (single scattering) contribution. Note the absence of cusp effects in 
the full result for E^ n . For Mf }l , there is no tb contribution to this order 
in the chiral expansion. Units are 10 _3 /M 7r +. 




Figure 6: Electric P-wave multipoles as a function of the pion excess en- 
ergy AW for photon virtuality k 2 = —0.1 GcV 2 . The solid (dashed) line 
shows the total (single scattering) contribution. Units are 10 _3 /M 7r +. 
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Figure 7: Longitudinal P-wave multipoles as a function of the pion 
excess energy AW for photon virtuality k 2 = -0.1 GcV 2 . The solid 
(dashed) line shows the total (single scattering) contribution. Units are 
10- 3 /M w +. 
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Figure 8: Magnetic P-wave multipoles as a function of the pion ex- 
cess energy AW for photon virtuality k 2 = —0.1 GcV 2 . The solid 
(dashed) line shows the total (single scattering) contribution. Units are 
lir z /M^. 
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w.f. 


NLO-500 


NLO-600 


NNLO*-500 


NNLO*-600 


«3 [GeV~ 4 ] 
< [GeV~ 4 ] 


4.010 
-5.925 


3.459 
-5.745 


4.832 
-5.895 


4.966 
-5.660 


a™ [GeV- 5 ] 


-34.49 


-36.45 


-29.86 


-27.51 



Table 2: Values of the fitted LECs for the various wave functions (w.f.). The 
values for 4 refer to the fits 2, whereas the corresponding a?? belongs to the 
respective fits 1. 



to the three-body terms. This was indeed done for the case of neutral pion photoproduction off deuterium in Ji6[ 
and we anticipate a similar result for the case under consideration (a complete fourth order calculation for the 
deuteron case can only be done when a similar investigation for single nucleon electroproduction is available and 
the already discussed inconsistencies have been resolved). We note in particular the cusp-like effects in certain 



P-wave multipoles due to the pion mass difference in the three-body contributions as discussed in Section 4.3 
A comparable multipole analysis of the data is not yet available. The multipoles extracted in Ref. || are based 
on the simplifying assumptions of constant S-waves and P-waves that solely depend on the pion center-of-mass 
momentum. A direct comparison of the multipoles obtained here with the ones of Q have therefore to be taken 
cum grano salis. Nonetheless we have performed the fits of type 1 by matching to the empirical value of \Ld\ 
to get a better handle on the theoretical uncertainties of our calculation. Note also that the recently proposed 
exact cancellation p5| between the single nucleon rescattering and the charge exchange three-body diagram 
at threshold is visible in the S-wave multipoles, the cusp effects from the ss and tb terms neatly cancel, cf. 
Fig. [|. However, this argument only affects a subset of graphs and does not lead to the conclusion that one is 
essentially sensitive to the single scattering amplitude, as reflected in our results. 

In Figs. [|,[l0]we show the differential cross sections for fits 1 and 2 employing the NNLO*-600 wave functions in 
comparison to the MAMI data H . These two lines corresponding to the two fit procedures can be considered as 
a measure of the theoretical uncertainty at this order. This uncertainty is comparable to the the experimental 
errors. The bell-shape behaviour of the differential cross sections at the higher values of the pion excess energy 
is similar to what is found in neutral pion photoproduction off protons and can be traced back to the large 
and delta-dominated third order P-wave LECs b Ptn , which are well described in terms of resonance saturation. 
Within large fluctuations, the data show more of a backward-forward angle asymmetry This feature might be 
better described when the P-waves have also been worked out to forth order, but it is fair to state that we do 
not find sizeable discrepancies between the data and the theoretical predictions. As in pion production off the 
proton, the S-wave multipoles are only dominant very close to threshold and the P-waves start to dominate at 
excess energies of a few MeV. 

The corresponding total cross sections as a function of the excess energy AW and of the photon polarization 
e are shown in Figs, [fl] and [l2] for fit 1 and in Figs. [ll] and [TJ for fit 2 and the NNLO*-600 wave function, 
respectively. We notice that for fit 1 with increasing excess energy and, in particular, with increasing photon 
polarization the data are systematically below the chiral prediction. Due to the fitting procedure, the slopes 
of the various curves for the Rosenbluth separation shown in Figs. [l2| are of course correct. On the other 
hand, the fitting procedure 2 gives an overall better description of the total cross sections, cf. Fig. |l3| with 
a somewhat too small longitudinal S-wave contribution, as most clearly seen in Fig. |l4], where again the 
Rosenbluth separation of the total cross section is plotted. These observation can further be sharpened by 
considering the dominant longitudinal multipoles as visualized in Fig. [IB], where the transverse and longitudinal 
threshold S-wave multipoles Ed and Ld are shown in comparison to the data for fits 1 and 2 using the NNLO*- 



600 wave function. We also note that Ed is slightly below the data from SAL 26 , whereas the prediction 
from |ll| was by the same amount above the data. This can be traced back to a variety of effects. First, we 
use slightly different input parameters (for the neutron) so that the single scattering contribution is somewhat 
reduced. Further, in contrast to Ref.|lj| we include the pion mass difference in the three-body contribution, 
which further reduces Ed by about 0.3 • 10 _3 /A J f Tr , compare Fig. || (the energy dependence of Eq ± is almost the 
same at the photon point fc 2 = 0). Also, our treatment of the Fermi motion (boost correction) is improved as 
compared to that paper and we thus have an additional reduction of Ed- Given that there are other fourth order 
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Figure 9: Differential cross section at AW — 0.5 MeV (left column) and AW = 1.5 MeV (right column) 
at three different values of the photon polarization for the NNLO*-600 wave function in comparison to 
the MAMI data §. Fit 1 (2): dashed (solid) lines. 
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Figure 10: Differential cross section at AW — 2.5 MeV (left column) and AW = 3.5 MeV (right column) 
at three different values of the photon polarization for the NNLO*-600 wave function in comparison to 
the MAMI data §. Fit 1 (2): dashed (solid) lines. 
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Figure 11: Total cross section as a function of AW for three different 
values of the photon polarization in comparison to the MAMI data [|| 
for fit 1 and the NNLO*-600 wave function. 




Figure 12: Total cross section as a function of the photon polarization 
e for four different values of the pion excess energy AW in comparison 
to the MAMI data 101 for fit 1 and the NNLO*-600 wave function. 
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Figure 13: Total cross section as a function of AW for three different 
values of the photon polarization in comparison to the MAMI data [|| 
for fit 2 and the NNLO*-600 wave function. 
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Figure 14: Total cross section as a function of the photon polarization 
e for four different values of the pion excess energy AW in comparison 
to the MAMI data |S| for fit 2 and the NNLO*-600 wave function. 
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effects, our result for the transverse threshold multipole is consistent with the data. In Table || we collect the 
S-wave cross section aod for the various wave functions and fit procedures. We remark again that the scaled S- 
wave cross section given in O was much too large because the dominant longitudinal multipole was not correctly 
represented. Thus, the dramatic difference between the CHPT prediction and the data has disappeared, and the 
overall description of the data is satisfactory but still needs to be improved. This will presumably be achieved 
when a complete fourth order calculation including the dominant isospin breaking effects has been performed. 
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Figure 15: Threshold multipoles \Ed\ and \Ld\ as a function of the pho- 
ton virtuality in comparison to the photon point data from SAL [ p6[ and 
the electroproduction data from MAMI |3| . The sign of the experimental 
result for Ld is taken to agree with the theoretical prediction. In our fits, 
the positive sign for Ld is preferred. Solid (dashed) lines: Fit 2 (1). The 
NNLO*-600 wave function is used. 



6 Summary and conclusions 

In this paper, we have studied neutral pion electroproduction off deuterium in the framework of chiral pertur- 
bation theory at and above threshold. The salient ingredients and results of this work can be summarized as 
follows: 

i) We have developed a general multipole decomposition for neutral pion production off spin-1 particles that 
is particularly suited for the threshold region and formulated in close analogy to the standard CGLN 
amplitudes for pion production off nucleons (spin-1/2 particles). Similar work was previously published 



ii) The interaction kernel and the wave functions are based consistently on chiral effective field theory. The 
kernel decomposes into a single scattering and a three-body contribution. We have chirally expanded the 
various contributions working to first non-trivial loop order 0(q 3 ), with the exception of the S-waves for 
the single scattering contribution. These have to be included to fourth order with one additional fifth 
order term [fill . 
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-k 2 [GeV 2 ] 


0.00 


0.02 


0.04 


0.06 


0.08 


0.10 


a od (NLO-500) [1] 


0.029 


0.031 


0.032 


0.035 


0.039 


0.046 


a od (NLO-500) [2] 


0.029 


0.022 


0.016 


0.014 


0.013 


0.013 


a od (NLO-600) [1] 


0.027 


0.029 


0.030 


0.032 


0.037 


0.043 


a od (NLO-600) [2] 


0.027 


0.020 


0.016 


0.013 


0.013 


0.014 


a od (NNLO*-500) [1] 


0.033 


0.035 


0.036 


0.039 


0.043 


0.048 


a od (NNLO*-500) [2] 


0.033 


0.024 


0.018 


0.014 


0.013 


0.012 


a od (NNLO*-600) [1] 


0.033 


0.035 


0.037 


0.039 


0.043 


0.048 


a od (NNLO*-600) [2] 


0.033 


0.025 


0.018 


0.015 


0.013 


0.013 



Table 3: S-wave cross section ao d in fib for the various wave functions 
(w.f.) and fit procedures [n] (n = 1,2 ) employed. 



iii) All parameters for pion production off the proton and the ones appearing in the three-body terms are 
fixed. The longitudinal neutron S-wave amplitude contains effectively two parameters, which we have 
determined by two different procedures. In the fits of type 1 we have fitted the fifth order parameter to 
the threshold multipole L d from Ref. || (and assuming resonance saturation to pin down the other LEC). 
The second procedure is based on a two parameter fit to the total cross section data from Ref. ||. All 
results are completely insensitive to the wave functions used, showing that this reaction is sensitive to the 
long-range pion exchange firmly rooted in the chiral symmetry of QCD. 

iv) The predicted differential cross sections are satisfactorily described for both fit procedures, although 
some systematic discrepancies for the higher values of the excess energy AW remain, see Figs. ^,[l0|. In- 
particular, for fit 1 the total cross section rises too steeply with pion excess energy. 

v) The calculated S- and P-wave multipoles exhibit a more complex pion energy and photon virtuality 
dependence as assumed in the fits of Ref. || . Within one standard deviation, the chiral predictions for 
the threshold multipoles \E d \ and \L d \ are consistent with the data at k 2 = p6| and k 2 = —0.1 GeV 2 ||. 

Clearly, the calculation presented here needs to be improved, in particular, the fourth order corrections to the 
P-waves and the three-body terms have to be included (note that similar work for the P-waves in neutral pion 
production off protons has only appeared recently Q). However, we have demonstrated that chiral perturbation 
theory can be used successfully to analyze pion electroproduction data off the deuteron which gives access to 
the elementary neutron amplitude. It would be very interesting to also have data at lower photon virtuality, 
which might also help to resolve the mystery surrouding the proton data at k 2 = —0.05 GeV 2 . 
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A Equivalence of multipole expansions 

In this appendix, we show the formal equivalence between the multipole expansion employed here and the earlier 
one developed in ||. Our transition amplitude has the form (we suppress here the phase factor related to the 
ip dependence) 

t m 'Xm{9)= (LnX + m-m'lm'\J\ + m)(lmL\\J\ + m)O^jY LntX+m _ m ,(9,cp = 0) (A.l) 

L„,L,J 
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in terms of the spherical harmonics. Arenhovel M works in the helicity basis and uses the rotation matrices 
<'<!.»■{<>■■ 

Cx, m {°)= E < L - 01 -m'\J-m'){l -mL\\J\-m)-^LO L L: X J d(_ m ^ m ,(9) . (A.2) 



We now show that 

1 

*m',A,m(0) = E *-M,A,-m(f) dl p ,-m'W ■ 

This can be proven simply by using the relation between the spherical harmonics and the d-functions, 



(A.3) 



L 



Y L ^x+rn-mi{0,<p = 0) = — = d^ +m _ m , „ (0) 



47T 



(A.4) 



and the following relation between rotation matrices 

E (?i ™i Ja '»2 | i m ) m , (0) m , (0) = £ (?i mi j 2 m 2 | j m' ) < m , (0) . (A.5) 



With this, the equivalence between Eqs. (A.l) and flA.2j) follows immediately. 



B Transformation matrices 



In this appendix, we collect the expressions of the various matrices appearing in Eqs. (2.13 2.14). For notational 
simplicity, we substitute the symbol for the pion angular momentum P^ by P in this appendix. Consider first 
Dl. The non-zero matrix elements are: 
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Dm 
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2^ 
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D 13 = 
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16 VL 

-* VTTT {p - 2+L - Pl) 

-8 + 16 P 
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L (P-i+l - Pi+l) L(P L - P 2+l ) = ~(LPl) ^ -(LP 1+l ) 
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-((3 + 5L + 2P 2 ) P_ 1+L + (1-2L) LP 1+L ) 
48v/ 1 £V^lT2l(l + 2L) 

VZVT+L {-P-i+l + Pi+l) 
> \ I -2/. " ' 

=i ((-3-2£ + 3L 2 + 2L 3 ) P_ 2+L + (3 + 4L - 6L 2 - 4L 3 ) P L + L (-2 + 3L + 2X 2 ) P 2+L ) 

v / Iv / TTTVl + 2L (-3 + 4X + 4L 2 ) 
^7TT2XPl ^ _ -g ((! + £) P-i+l + LP 1+l ) 

VlVTTl ' 74 v / Lv / TTTVl + 2L 

-((-3 + P + Z 2 ) (P i+l - Pi+l)) 
24 VIV1TIVTT2T 

-((3 + 2L) 2 (-1 + i 2 ) P-2+l) + (-9- 14L + 12L 2 + 8L 3 ) P L + (l-2£) 2 £ (2 + L) P 2+L 

48VI\/TTIV1 + 2L (-3 + 4X + 4L 2 ) 
VIVT+T (-P-i+l + Pi+l) d - (VTT2I Pl) 

24V1 + 2L ' 78 16 VI VTTI ' 

-((3 + 5L + 2L 2 ) P_ 1+L + (1-2L) L P l+L ) 

48 VI vTTX V1T2T 

yiyr+Tv3 + 2L (p_ 1+L -Pi+ £ ) 
8 + 16L 
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D 82 = 



D 83 = 



D 8e = 



D 87 = 



1+L 



((-3 -2L + 3L 2 + 2L 3 ) P_ 2+L + (3 + AL - 6L 2 - 4L 3 ) P L + L (-2 + 3L + 2L 2 ) P 2+L ) 



V3 + 2T (-1 + 4L 2 ) 



L 



1 + L 



V3 + 2LPl , D 



84 



£ 
1+L 



V3 + 2L ((1 + L) P- 1+L + L P 1+L ) 



16 L 



(-3 + L + L 2 ) (P_ 1+L -P 1+L ) 



24 (1 + 2L) V3 + 2L 



^ (3 + 2L) 2 (-1 + i 2 ) P- 2+L + (9 + 14L-12L 2 -8L 3 ) P L -(1-2L) 2 L (2 + L) P 2+L 



D. 



89 



A)5 = 



A>7 = 



P 9 9 



48(3 + 2L)^ (-1 + 4L 2 ) 



X^TTTL (-P_ 1+L +P 1+£ ) 
24 (1 + 2L) V3 + 2L 



P 
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1+L 



Pl 



16V3 + 2L ' 



^ ((3 + 5L + 2I 2 ) P_ 1+L + (1-2L) iP 1+L 



48 (1 + 2X) V3 + 2L 



LVT+LV3 + L (P-i+l -Pi+l) 
8 (1 + 2L) V3 + 2L 



Pfl6 = 



(1 + L)^V3+X(Pl-P2+l) 



8(3 + 2i) 



VT+T (2 + L) V3+T (Pi +L - P 3+£ ) 
8 V3 + 2 Z (5 + 2L) 



A>8 = 



VT+lV3 + lp l 



8V3 + 2X 



Vi + L V3 + L Pi+l 



8V3 + 2X 



(B.l) 

Here, the Pl are the conventional Legendre polynomials that depend on z = cos 9. Note that L is positive 
definite so that Legendre polynomials with negative index have to be understood as zero. The matrix El has 
no zero elements and takes the form 



( 



\ 



^■V-i+lp l 

V2 

-y x/2+TP L 
V2 



v^T+TP-i+z, x'-l+L (l+L) (P_ 1+I ,-P 1+I ,) x^T+Tq+L) (f-^L-Ptj 



72 

| V^+£P-n 
V2 

I VT p 1+i , 

V2 



\ 



(-!+£) % /I+£(P-i + i,-Pi + L ) Vl+T(P- 2+I ,-P I ,) 

8V2(1+2L) 8v / 2(-l+2£) 

-(■/£ (2+£) (P-i + L -Pi +I ,)) -(yT(2+£) (Pj,-P 2+ i.)) 

8v^(l+2£) 8V2(3+2£) 

-(£ V2+L(P_ 1+L -P 1+L )) -(LV2+L(P L -P 2+L )) 

8V2(l+2£) 8\/2(3+2£) 



(B.2) 



We now turn to the matrices Gl„ and Hl^ appearing in the inverse transformation, Eq. (2.14). It is most 
convenient to express these with the help of Clebsch-Gordan coefficients. For that, we employ the P-symbols 



P 



L,L W ,J 
X,m' ,m 



= (Lit A + m — m 1 m \ J A + m)(l mL A| J A + m) 



(B.3) 



as they appear also in the multipole expansion of the T-matrix. We define the transverse and the longitudinal 
vectors T(L, J) and L(L, L % , J), respectively, in terms of their components: 



T(L, L„, J)i 

T(L, L„, J)2 
T[L, L„, J)3 



4VT+2LVT+2L~ [ D^'-i + Di'oo' J + ^ii 



i£,£,r,J 



P. 



(1) 



3V1 + 2J A /L w (l + L ff ) 
2iV2 + 4LVl + 2i ff (p^lj + P 1 J ; L 1 7r ' J ) P, 



(2) 
£„ 



VT+2J^/L V (1 + L„) y/-2 + L„ + Ll 
J V2 + 4lVl + 2^(Pf; L r 1 7 + D^> J )P L „ 



V1 + 2J 
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+ 



T(L, L v , J)i 



2iy/l + 2LV1 + 2L 7I z(D^ I :^ 1 - Dffi' J )p£} 

VT+2Jy/L v (l+L„) 

»V2 + 4iyi + gig-i + z 2 )(pf;fr; J + rffirM? 

VT+2JV L A1 + L„)y/-2 + L^ + Lk 2 
2^^/l + 2LV1 + 2L^{D^l'^ - D^J' J )P { ^ 
VT+2J^/L^(r+L^ ' 

2 VI + 2Ly/l + 2Ln / n L : L x ,J ort L,L„,J n L,L„,J , n L,L x ,./^ 



T(L,L n ,J) 7 
T(L, L n , J)s 



L(L, L n , J) 2 
Ltt, J)3 



+ 



+ 



T(L,L V ,J) 6 — 



+ 



+ 



T(L,L v ,J)g — 



L(L,L v ,J)i — 



+ 



(2) 



2Vl + 2LVl + 2i w (l + 3z 2 )L>[ '^I'lPjP 

VI + 2Jy/L v (l + i 7 r)v / -6 + L T + ZorV-2 + ivr + ^ 

4V2 + 4LV1 + 2L„(£>ffi/ - ggr^Pff 
VI + 2Jy/L„(l + L v )y/-2 + L v + L 2 

16V1 + 2LV1 + 2L lx zD^l{P [ l ) 



VI + 2Jy/L„(l + ^tt) V-6 + L n + L^ 2 V 2 + L ff + L^ 2 
8V1 + 2LV1 + 2L v D^fp^ 



VI + 2JV^7r(l + L„)y/-6 + L n + L 2 yJ-2 + L n + L„ 
2V2 + 4IV1 + 2^(^/0^/ - d[£« J )Pl„ 
V1 + 2J 

4VTT2LV1 + 2L^zDi,t-i P L« 



VTT2JV^( 1 + ^) 



2V2 + 4lVT+2i;(-l + z 2 )(£>f/iT,b J - I^rM* } 



VI + 2Jy/L„(l + L v )y/-2 + i w + I^ 2 

4VTT2IV1 + 2L„z(-i + z 2 )^f;fr;/pff 

VI + 2Jy/L*(\ + irr) V"6 + L w + L„ 2 y/-2 + L^ + L 2 
-4VTT2lyi + 2UB L 1 ^j(P { l] 

VT+Ij^l^i + Lk) 



4VTT2IVT+2I7(-1 + z 2 )D^fP^ 



VI + 2Jy/L*(l + L*)y/-6 + Lit + L 2 ^/-2 + L n + L 2 
2iV2 + 4lVl + 2L v D^T J Pl„ 
V1 + 2J 

2 ? VrT2lyi + 2L v z(D^l- J - D^ ,J )Pl} 

VT+27V L -( 1 + i -) 

-2iyT+2Wl + 2L^D^; L J' J - D^ J )P^ 

VTT2jyi7TTX7) 

-4VTT2iyT+2i7(j L b L r J + Dt L r J ) p ^ 

VT+2Jy/L„(i + L v ) 
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V2 + 4LV1 + 2L 7V zD^^P^ 2) 



VI + 2Jy / L 7T (l + L n )y/-2 + L n + L 

\/l + 2J V /L 7r (l + L 7r )v / -2 + Z/Tj- + L^ 2 
:ivative of the 

expressed in terms of the following vectors 



-4V2 + + 2L 7t D^^P {2) 

L(L,L K ,J) 4 = * L ", , ■ (B.5) 

where P^\z) is the n th derivative of the Legendre polynom Pl(z). The matrices Gl„ and if^ can then be 



with 



Gl„ = (G l ,...,G 9 ) , (B.6) 
H Lk = (#!,..., if 4 ) , (B.7) 



Gi = T(i 7r — 1, L,,., L,r — 1) , Gi = T{L-n — 1, L^, L x ) , G3 = T(i 7r + 1, Lj,-, L„.) , 

G4 = T(i 7r + 1, L,,-, L w + 1) , G5 = T(L W — 2, Ljr, — 1) , Gg = T(L V , — 1) , 

G7 = T(L n , Ljr, L„.) , Gg, = T{L 7! ,L 7T ,L 7T + 1) , Gg = T^L^ + 2, L ff , + 1) , (B. 



and 



J?i — L(Ln — 1, L^, — 1) , — L(L 7r — 1, Lj,-, Ltt) , 

) , if 4 = + 1, + 1) . (B.9) 

C Two-body to three-body center-of-mass 

In this appendix we sketch the derivation of the transformation from the 7-d center-of-mass (COM) system to 
the 7- TV COM, extending the considerations given in We are interested in the kinematics of the process 
j*NiN2 — > ttNxN2, where the nucleons, N% and N2, are sewn to the deuteron wavefunctions. Our 3-body 
corrections are evaluated in the 7-d COM whereas the single scattering corrections which take into account 
the scattering of the photon on the individual nucleons have been calculated in the j-N COM. It is therefore 
necessary to construct the Lorentz transformation which boosts the single-scattering corrections to the j-d 
COM. 

Let p be some four-vector in the 7-d COM and p* the corresponding four-vector in the COM of the (second) 
single nucleon. These are related by the Lorentz transformation p* — A(u)p with u — ue x the velocity. The 
vector p transforms as 

Pa \ _ I 7 ~h 
P*\ J \ -Pi 7 
and the transverse directions are of course unaffected. This gives 




(C.l) 



£=^L±*L = -_£—, (c.2) 

k + P2Q k + P20 



so that e x = —pi . The photon energy-momentum four-vector thus transforms as 

fc„ = 7 (fc + file ■ pi ) , 



k* = k- (fc.fc(l-7)-7/3fcoj pi . (C.3) 
Expanded in powers of 1/m, this reads 

fco = fco + -fc-Pi+C(l/m 2 ) , 
m 

k* = k+-k Q pi + 0(l/m 2 ) , (C.4) 
m 
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and similarly for the pion energy and three-momentum (go = ^,<f)- One also has to transform the photon 
polarization vector. This is most easily done if one uses the following gauge-invariant form of the 7* N — > ir°N 
transition amplitude, 

6 

M„=x;o»,iic,i > (c.5) 

i=l 

with 

= e* ■ (fc* x <f) , SS , 2 - £* ■ (fc* xf)S- (k* x <f ) , O ss , 3 = t ■ (k* x (k* x £)) , 

(C.6) 

Since we work in the Coulomb gauge in the 7-d COM, the full Lorentz transformation for the polarization vector 
is given by 

£o = 7/3e-pi, e* = e - e -pi (1 - -f)pi . (C.7) 
Expanded in powers of 1/m, one observes that only the time-component is modified to leading order, 

e* = — e-px +C(l/m 2 ) , e * = e + C(l/m 2 ) . (C.8) 
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